Optimal Lβ\mathfrak{L}^{\beta}-Control for the Global Cauchy Problem of the Relativistic Vlasov-Poisson System

Recently, M.K.-H. Kiessling and A.S. Tahvildar-Zadeh proved that a unique global classical solution to the relativistic Vlasov-Poisson system exists whenever the positive, integrable initial datum is spherically symmetric, compactly supported in momentum space, vanishes on characteristics with vanishing angular momentum, and for $\beta \ge 3/2$ has $\mathfrak{L}^{\beta}$-norm strictly below a positive, critical value $\mathcal{C}_{\beta}$. Everything else being equal, data leading to finite time blow-up can be found with $\mathfrak{L}^{\beta}$-norm surpassing $\mathcal{C}_{\beta}$ for any $\beta >1$, with $\mathcal{C}_{\beta}>0$ if and only if $\beta\geq 3/2$. In their paper, the critical value for $\beta = {3}/{2}$ is calculated explicitly while the value for all other $\beta$ is merely characterized as the infimum of a functional over an appropriate function space. In this work, the existence of minimizers is established, and the exact expression of $\mathcal{C}_{\beta}$ is calculated in terms of the famous Lane-Emden functions. Numerical computations of the $\mathcal{C}_{\beta}$ are presented along with some elementary asymptotics near the critical exponent ${3}/{2}$.Comment: 24 pages, 2 figures Refereed and accepted for publication in Transport Theory and Statistical Physic

Asymptotically complexity diminishing schemes (ACDS) for kinetic equations in the diffusive scaling

In this work, we develop a new class of numerical schemes for collisional kinetic equations in the diffusive regime. The first step consists in reformulating the problem by decomposing the solution in the time evolution of an equilibrium state plus a perturbation. Then, the scheme combines a Monte Carlo solver for the perturbation with an Eulerian method for the equilibrium part, and is designed in such a way to be uniformly stable with respect to the diffusive scaling and to be consistent with the asymptotic diffusion equation. Moreover, since particles are only used to describe the perturbation part of the solution, the scheme becomes computationally less expensive – and is thus an asymptotically complexity diminishing scheme (ACDS) – as the solution approaches the equilibrium state due to the fact that the number of particles diminishes accordingly. This contrasts with standard methods for kinetic equations where the computational cost increases (or at least does not decrease) with the number of interactions. At the same time, the statistical error due to the Monte Carlo part of the solution decreases as the system approaches the equilibrium state: the method automatically degenerates to a solution of the macroscopic diffusion equation in the limit of infinite number of interactions. After a detailed description of the method, we perform several numerical tests and compare this new approach with classical numerical methods on various problems up to the full three dimensional case

Wavelet approximations of a collision operator in kinetic theory

International audienceApproximations par ondelettes d'un opérateur de collision en théorie cinétique. Cette Note présente un schéma conservatif et entropique basé sur l'utilisation des ondelettes pour l'équation de Fokker-Planck-Landau isotrope qui modélise l'évolution des particules chargées dans un plasma. La présente approche possède à la fois les propriétés des schémas aux différences finies (conservation et entropie) et celles des méthodes spectrales (essentiellement la précision) qui sont développés dans la littérature. De plus, l'approche ondelette fournit un algorithme rapide d'évaluation de l'opérateur de collision. Ce travail constitue un premier pas vers le développement d'approximations par ondelettes d'opérateurs de collision plus complexes

A new variational approach to the stability of gravitational systems

We consider the three dimensional gravitational Vlasov Poisson system which describes the mechanical state of a stellar system subject to its own gravity. A well-known conjecture in astrophysics is that the steady state solutions which are nonincreasing functions of their microscopic energy are nonlinearly stable by the flow. This was proved at the linear level by several authors based on the pioneering work by Antonov in 1961. Since then, standard variational techniques based on concentration compactness methods as introduced by P.-L. Lions in 1983 have led to the nonlinear stability of subclasses of stationary solutions of ground state type. In this paper, inspired by pioneering works from the physics litterature (Lynden-Bell 94, Wiechen-Ziegler-Schindler MNRAS 88, Aly MNRAS 89), we use the monotonicity of the Hamiltonian under generalized symmetric rearrangement transformations to prove that non increasing steady solutions are local minimizer of the Hamiltonian under equimeasurable constraints, and extract compactness from suitable minimizing sequences. This implies the nonlinear stability of nonincreasing anisotropic steady states under radially symmetric perturbations

Asymptotic Preserving scheme for a kinetic model describing incompressible fluids

International audienceThe kinetic theory of fluid turbulence modeling developed by Degond and Lemou in [6] is considered for further study, analysis and simulation. Starting with the Boltzmann like equation representation for turbulence modeling, a relaxation type collision term is introduced for isotropic turbulence. In order to describe some important turbulence phe-nomenology, the relaxation time incorporates a dependency on the turbulent microscopic energy and this makes difficult the construction of efficient numerical methods. To investi-gate this problem, we focus here on a multi-dimensional prototype model and first propose an appropriate change of frame that makes the numerical study simpler. Then, a numerical strategy to tackle the stiff relaxation source term is introduced in the spirit of Asymptotic Preserving Schemes. Numerical tests are performed in a one-dimensional framework on the basis of the developed strategy to confirm its efficiency

Orbital stability of spherical galactic models

International audienceWe consider the three dimensional gravitational Vlasov Poisson system which is a canonical model in astrophysics to describe the dynamics of galactic clusters. A well known conjecture is the stability of spherical models which are nonincreasing radially symmetric steady states solutions. This conjecture was proved at the linear level by several authors in the continuation of the breakthrough work by Antonov in 1961. In a previous work, we derived the stability of anisotropic models under {\it spherically symmetric perturbations} using fundamental monotonicity properties of the Hamiltonian under suitable generalized symmetric rearrangements first observed in the physics litterature. In this work, we show how this approach combined with a {\it new generalized} Antonov type coercivity property implies the orbital stability of spherical models under general perturbations

Stable ground states for the relativistic gravitational Vlasov-Poisson system

We consider the three dimensional gravitational Vlasov-Poisson (GVP) system in both classical and relativistic cases. The classical problem is subcritical in the natural energy space and the stability of a large class of ground states has been derived by various authors. The relativistic problem is critical and displays finite time blow up solutions. Using standard concentration compactness techniques, we however show that the breaking of the scaling symmetry allows the existence of stable relativistic ground states. A new feature in our analysis which applies both to the classical and relativistic problem is that the orbital stability of the ground states does not rely as usual on an argument of uniqueness of suitable minimizers --which is mostly unknown-- but on strong rigidity properties of the transport flow, and this extends the class of minimizers for which orbital stability is now proved

An asymptotic preserving scheme for the Kac model of the Boltzmann equation in the diffusion limit

International audienceIn this paper we propose a numerical scheme to solve the Kac model of the Boltzmann equation for multiscale rarefied gas dynamics. This scheme is uniformly stable with respect to the Knudsen number, consistent with the fluid-diffusion limit for small Knudsen numbers, and with the Kac equation in the kinetic regime. Our approach is based on the micro-macro decomposition which leads to an equivalent formulation of the Kac model that couples a kinetic equation with macroscopic ones. This method is validated with various test cases and compared to other standard methods

A non-variational approach to nonlinear stability in stellar dynamics applied to the King model

In previous work by Y. Guo and G. Rein, nonlinear stability of equilibria in stellar dynamics, i.e., of steady states of the Vlasov-Poisson system, was accessed by variational techniques. Here we propose a different, non-variational technique and use it to prove nonlinear stability of the King model against a class of spherically symmetric, dynamically accessible perturbations. This model is very important in astrophysics and was out of reach of the previous techniques